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Abstract 

The proximity of Z + (4430) to the D*Di threshold suggests that it may be a D*Di molecular 
state. The D*Di system has been studied dynamically from quark model, and state mixing effect is 
taken into account by solving the multichannel Schrodinger equation numerically. We suggest the 
most favorable quantum number is J p = CP, if future experiments confirm Z + (4430) as a loosely 
bound molecule state. More precise measurements of Z + (4430) mass and width, partial wave 
analysis are helpful to understand its structure. The analogous heavy flavor mesons and Z^ + 
are studied as well, and the masses predicted in our model are in agreement with the predictions 
from potential model and QCD sum rule. We further apply our model to the DD* and DD* system. 
We find the exotic DD* bound molecule doesn't exist, while the 1 ++ DD* bound state solution can 
be found only if the screening mass \i is smaller than 0.17 GeV. The state mixing effect between the 
molecular state and the conventional charmonium should be considered to understand the nature 
of X(3872). 

PACS numbers: 12.39.Jh, 12.40.Yx,13.75.Lb 
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I. INTRODUCTION 



In the past years many new mesons have been discovered through B meson decays. 
Recently the Belle Collaboration has reported a narrow peak in the 7r + ip' invariant mass 



spectrum in B — > K^ip' with statistical significance greater than 7a This structure 
is denoted as Z + (4430). The Breit Wigner fit for this resonance yields the peak mass 
M = 4433 ± A(stat) ± l(syst) MeV and the width T = Ut\l(stat)^n(syst) MeV. The 
product branching fraction is determined to be £>(B — > KZ + (4430)) -£>(Z + (4430) — > ir + ifj') = 
(4.1±1.0(stat)±1.3(syst)) x 10" 5 . Since the G-parity of both vr+ and if/ is negative, Z+(4430) 
is a isovector with positive G-parity. However, Z + (4430) is far from being established, no 
significant evidence for Z~(4430) has been observed neither in the total J/tpir" or ip(2S)7r~ 
mass distribution by the Babar Collaboration 

There are already many theoretical investigations for the possible structures and the 



properties of Z + (4430) 
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Because it is 



very close to the threshold of D*Di(2420), and the width of Z + (4430) is approximately the 
same as that of D 1 (2420), it is suggested that Z+(4430) could be a D*Di(2420) molecular 



state 



B,y,y, 



ii 



1 61 ] . Other interpretations such as tetraquark state [5|, |8j or a cusp in the 
D*Di channel [l4] are proposed as well. In Ref. jjj, we suggested how to distinguish the 
molecule and the tetraquark hypothesis, and Z + (4430) as a D*D! molecule was studied from 
the effective field theory. In Ref. fill . the authors investigated dynamically whether 
Z + (4430) could be a S-wave D*Di or D*!)^ molecular state by one-pion exchange and a 
exchange. 

In principle, nothing in QCD prevents the formation of nuclear-like bound states of 
mesons and speculation on the existence of such states dates back thirty years 19| . Tornqvist 
suggested that two open flavor heavy mesons can form deuteron-like states due to the strong 



it exchange interaction 



20| | . and the monopole form factor is introduced to regularize the 



interaction potential at short distance. In Ref. 21], the author investigated the possible 
heavy flavor molecules base on long distance one pion exchange and short distance quark 
interchange model. However, the dynamics of hadronic molecule is still unclear so far. In 
this work, we will dynamically study Z + (4430) and analogous heavy flavor states Z^ and 
Zj. + from quark model. We shall discuss the interaction between two hadrons at the quark 
level instead of at the hadron level. The effective interactions between quarks including the 
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screened color-Coulomb, screened linear confinement and spin-spin interactions are employed 
to describe the interactions between the components of the interacting hadrons. 

This paper is organized as follows. In section II, the canonical coordinate system and 
the effective interactions are introduced. We give the details of the evaluation of the matrix 
element in section III. In section IV, the D*T)i system coupled with D*D 2 is studied, and 
the possible structure of Z + (4430) is discussed. In section V, the analogous heavy flavor 
states Z^ b and Z^, + , DD* and DD* systems are investigated, the static properties such as the 
mass and the root of mean square radius etc. are calculated. We present our conclusions 
and some relevant discussions in section VI. Finally the spatial matrix elements involved are 
given in the Appendix. 

II. CANONICAL COORDINATE SYSTEM AND THE EFFECTIVE INTERAC- 
TIONS 




FIG. 1: Canonical coordinate system for the four quark system, where black circle denotes quark 
and empty circle denotes antiquark. 

The coordinate shown in Fig. [I] is taken as the canonical coordinate system, which defines 
the asymptotic states. The relevant coordinates of this system can be expressed in terms of 
i"i3, i" 24 and r as follows, 

m 3 7714 
T\2 = fi3 — r?4 — r 

mi + r«3 7712 + 777.4 
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m 3 m 2 

~ r 13 H ; r 24 — r 



r u = 

mi + m 3 m 2 + m 4 
mi m 4 

1*32 = I"l3 — 1*24 — T 

mi + m 3 m 2 + m 4 
mi m 2 
T34 = ; r i3 H ; r 24 — r 



(1) 



mi + 777.3 ^2 + 

where m 1; m 2 , m 3 and m 4 are respectively the masses of constituents 1, 2, 3 and 4. The 
relative position r between the center of mass of the two mesons is 

m 1 m 2 r 12 + mim 4 ri 4 + m 2 m 3 r 32 + m 3 m 4 r 34 

r = 7 w s ( 2 ) 

(mi + m 3 )(m 2 + m 4 ) 

As is shown in Eq.flT]), we can compactly represent the coordinate r»j in terms of 1*13, r 24 
and r as follows 

r« = /a(^>i3 + /b(u>24 - r, z G A and j E B (3) 
The parameters and are listed in Table [B 
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TABLE I: The values for the parameters /a(*j) an d /b(*j)- 

In the above canonical coordinate, the Hamiltonian for this system, including the relative 
motion and the interaction between two mesons, is split into 

1 



H = H (A(13)) + H (B(24)) 



2/iAB 



-Vi + Vr 



(4) 



where jTo(^(13)) and Hq(B (24)) are respectively the Hamiltonian for the two mesons A and 
B, which contains the kinetic term and all the interactions within each meson. \xab is the 
reduced mass [Iab — Ma+Mb ' ^ ne ^ mr< ^ term — 2 ^ b is the kinetic energy operator of the 
relative motion. The interaction potential Vi is the sum of two-body interaction between 
quarks in the mesons A and B, 



(5) 
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The phenomenological interaction between a quark and an antiquark in a single meson ( e.g. 
A and B ) is reasonably well known, it is described by the short distance one- gluo n exchange 



interaction and the long distance phenomenological confinement interaction 

A(z') ra, 
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(6) 



where a s is the strong coupling constant, b is the string tension, and rrij are the masses 
of the interacting constituents. For an antiquark, the generator A/2 is replaced by -A T /2. 

Since we mainly concentrate on the molecular states comprising two heavy flavor mesons 
in this work, and the molecule is generally weakly bound. Therefore the separation between 
the two mesons in the molecule is rather larger than the average radius of the individual 
meson, and the two mesons interact mainly through two gluons exchange processes 



24 



25|], 



which results in the color van der Waals interaction. By compa ring with the van der Waals 
interaction between electric dipoles in QED, the author in Ref. 26| introduced the effective 
charges for quarks and antiquarks to describe the color van der Waals interaction between 



two mesons. The effective charges for quark and antiquark respectively are C„ 



2N C 



and Cn 



N?-l 



y 2N , here N c is the number of color with N c = 3 in QCD. It is remarkable 
that the effective charge correctly describes the interaction between quark and antiquark in 
an individual meson as well. The effective charge is also consistent with the Lattice QCD 
results, which found the nonperturbative potential between a quark and an antiquark in 
different representations is proportional to the eigenvalue of the quadratic Casimir operator 



27]. 



Different from the interactions between quarks in a single meson, as the interaction 
between the constituents in a molecule takes place at large distances, we are well advised 
to use a screened potential to represent the effects of dynamical quarks and gluon 28[. A 
simple way to incorporate the screening effect is to replace k 2 in the Fourier transformation 



of the interaction potential by k 2 + fi 2 



26 



29], where /i is the screening mass. With 



the effective charge and the screening effect in mind, in momentum space, the effective 



interaction potential (k) between two quarks in the mesons A and B is 



V, e / f (k) 



i.i 



Ana* 
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k 2 + /i 2 (k 2 + /i 
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(7) 



The effective interaction in coordinate space V£ (r 



2>mirrij 

is the Fourier transformation of 
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(k) 



(8) 



Therefore in coordinate space the effective interaction (fy) is 



V; 



eff. 



'J 
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(9) 
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In this work, we will use the above effective interaction Vf/H^ij) to study the possible 
heavy flavor molecules dynamically. Comparing with Ref. 26[, we have introduced the 
spin-spin interaction in addition to the screened color-Coulomb and the screened linear 
confinement interactions. In the light quark hadrons, the spin-spin hyperfine interaction 
makes the dominant contribution to the hadron-hadron interactions 30j. Whereas, for the 



heavy flavor mesons, the hyperfine interaction contribution is smaller due to the large heavy 
quark mass 3lJ, |3_2j. Therefore we expect that the contribution of spin-spin interaction 
should be smaller than those of the screened color-Coulomb and screened linear confinement 
interactions. However, the spin-spin hyperfine interaction may play an important role when 
we study the dynamics of molecular state, since the binding energy of molecular state is 
usually rather small. On the other hand, if we neglect the spin-spin interaction, the static 
properties of the molecule, such as the binding energy and the root of mean square radius 
(rms) etc., would be independent of the spin of the molecular state, which contradict with 
the experimental observations for deuteron. 

As a result of the residual interaction Vi between two mesons, at short distance the 
mesons may excite as they interact, and they could be virtually whatever the dynamics 
requires. This means that we need to consider the state mixing effect. It has been shown 
that the state mixing effect plays an important role in obtaining the phenomenologically 
required potential, when we study the nucleon-nucleon and nucleon-antinucleon interactions 



from the chiral soliton model 
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34j. The eigenvalue equation for the system is 



(H-E)\V) =0 



(10) 



where E and are respectively the eigenvalue and the corresponding eigenfunction. If 
there is no residual interactions Vj between A and B, the eigenfunction of the total system 
would simply be the product of A meson's wavefunction and B meson's. Consequently it is 
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natural to expand the eigenfunction in terms of the model wavefunctions 

\y) = Y,H*U®AAB)) (n) 

a' 

where ip(r) a is the relative wavefunction between the mesons A and B, and \<& a (A, B)) = 
\$a)\®b) denotes the intrinsic state of the two mesons, which will be mixed under the interac- 
tion Vj. The wavefunction |$^) satisfies the Schrodinger equation (H (A(13)) — Ma)\$a) = 
0, where \<&a) depends on the relative coordinate r 13 , and similarly for Inserting wave- 

function |\&) into the eigenequation Eq. fllOp . multiplying by (<& a | and integrating over the 
internal coordinates, we obtain 

V* + V Iaa (r) + E a - E)^ a (v) = -J2 Vw(r)^ a <(r) (12) 



Where E a = Ma + Mb is the energy eigenvalue of channel a. Vj aQ /(r) = ($ a \Vi\§ a i) is the 
matrix element of the interaction potential V/, it is a function of the relative coordinate r, 
the intrinsic coordinates r 13 and r 2 4 have been integrated out. There is clearly one equation 
for each state a, and they are coupled each other by the terms on the right-hand side. 
It is important to notice that all the transitions represented by the right hand of Eq. (112j) 
contribute coherently. If \E a — E a t\ » \Vi aa /(r)\ with a ^ a', then the coupled channel 
Schrodinger equation Eq. (fl2|) is reduced to the single channel Schrodinger equation 

( " + V; aa (r) + E a - E)Mr) = (13) 

where Vj aa {r) is the effective interaction potential 

VUr) = V Iaa (r) ~ E (14) 
Eq. lflBl) and Eq. (fl4|) are exactly the results of the second order perturbation theory to deal 



33 



341 . However, if 



with the state mixing effect, and this simplification is widely used [2i 
\E a — E a i\ is rather small or of the same order comparing with |Vr aa '(r)|, we have to solve 
the coupled channel Schrodinger equation exactly. Although in principle we should solve 
the infinite set of equations implied by Eq. (1121) . in practice we only need to concentrate on 
the nearly degenerate channels, which is a good approximation. 
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III. EVALUATION OF THE MATRIX ELEMENT Vw(r) 



For a system consisting of two mesons A and B with total angular momentum J and the 
third component J z , its wavefunction is written as 



\^(A,B))=J2{[(XAXBf(MB 
S.L 



[(Xa^aY a (xb^bY b ] J ' Jz ) [{xaXb)\M 



^ SLJ^Jfi 

S,Sz ,L 



Sa Sb S 
La Lb L 
Ja J_b J 



) (S, S 2 ; L, J, - S Z |J, J z )\(xaXb) S ' Sz )\(Mb) 



\ L, J z — S2 



(15) 



where S = V2S+T, x is the spin wavefunction, and tp represents the spatial wavefunction. 
Sa, La and denote respectively the spin, the orbital angular momentum and the total 
angular momentum of meson A with similar notations for the meson B. From Eq.([5l) and 
Eq.Q, it is obvious that each term of Vj can be factorized into the spatial and spin relevant 
part, consequently the interaction potential Vi can be re- written as 

3 



vi= E E^K (fe) (%)K (fe) 

iSAj'S-B k=l 



(16) 



where the superscript (k) represents respectively the screened color Coulomb, screened linear, 
and spin-spin interactions for k =1, 2, 3. Concretely, V} 1 ' = V} 2 ^ = 1, V} 3 ^ = Sj • Sj, and 
the spatial part V} k \rij) can be read from Eq.([9]) straightforwardly. Therefore the matrix 
element Vi aa >(r) is the sum of twelve terms, and each term is of the form 
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Sb 
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i',S' z ,L' 




















Ja 


Jb 


J 

) 




J A' 


Jb' 


J' 



x (s, s 2 ; l, h - s 2 | j, j 2 )(s', s; ; u, j; - s'jj', j' g ){{M B >f' J *- s W fe) l(^B) L ' J *~ s *) 



x((XA>XB>f> S '>\V^\(XAXB) S > S >) 



(17) 



It is obvious that both the spatial matrix element ((4>A'4>B') L '' J ' z ~ S ' z \Vr k ^\(' l l J A'4>B) L '' !z ~ Sz ) an d 
the spin matrix element ((xa'XbO S '' S HK^I(XaXb) S ' Sj! ) are needed. Firstly we consider the 
spatial matrix element 

((^ fl ')^-^|vW(ry)|(^) L ' J '- s '> = E (l^,l^;Lb,,Lb^|l', j; - s;> 

LAz,Lbz,L a i z ,L b i z 

(L a ,L A2 ;Lb,Lb 2 |L, J z - S z )<^'' L ^(r 13 )^^ (U 

8 



where 



(^A^'' LA ' z ( r i3)i ) B^'' hB ' z (' r 24 : )\V} k \r ij)|^ j4,L ' 4z (ri 3 )^ s ' Lsz (r24)) 
= (L A ',L A i z ; L B >, L B i z \V} k \rij)\L Al L Az ; L B , L Bz ) 



d 3 r 



d 3 r 24 U L A t" LA '*(r 13 



B l r 24y 



xV} k \f A (ij)r 13 + f B (ij)r 2i -r) 



(19) 



In this work, the spatial wavefunctions are taken as the simple harmonic oscillator wavefunc- 
tions, which is a widely used approximation in the quark model calculations. The integral 



in 



q. (TT9l) can be evaluated analytically in coordinate space following the procedures in Ref. 



35( | . On the other hand this integration can be performed in momentum space as well, then 



the calculation will be greatly simplified 
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32] 



(W, La'z] W, L B ' Z 

d 3 p 



/dp 
e~ tPT PL A ,,L A , z ;L A ,h A AfA^j)p} PL Bl ,L B , z ;L B ,L B AfB{l])p}V {k) {p) 



(20) 



where 



d 3 r 13 e 4p - ri3 



PL A ,,L A , Z ;L A ,L A AP) ~ 

^)(p) = /A jJ -e— «^)(r,) 



V^' L ^(r 13 ) 



(21) 



Note that V^(p) can be read from Eq. (l7|) directly. For the given quantum numbers L A , L Az 
etc, the above integral can be straightforwardly calculated although it is somewhat lengthy, 
and the matrix elements involved in our calculation are listed in the Appendix. 

Next we turn to the spin matrix element ((xa'Xb') S '' S ' z IK^' 1 \(XaXb) S,Sz ) • We denote the 
spin of the constituents 1, 2, 3 and 4 by Si, s 2 , s 3 and s 4 respectively. In the present work, the 
constituent is quark or antiquark, consequently we have si = s 2 = s 3 = s 4 = \. We would 
like to recouple the constituents so that the spin operator Vj k '(k =1,2,3) matrix elements 
can be easily calculated. We have 



KxaXb) S ' Sz ) = |[(s 1 s 3 )S a (s 2 s 4 )S b ] S ' S2 ) 

si s 3 S A 

— X] Si 2 s 34 s A Ss 

Sl2,S34 



s 2 s 4 S B 
Si 2 S 34 S 



M[(sis 2 )S 12 (s 3 s 4 )S 34 ] s ^) 
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E (-1) Sb - S2 ~ S4 Si 4 S 32 SaSb 

Sl4,S32 



si s 3 S A 
s 4 s 2 S B 

Sl4 S32 S 



M[(sis 4 )Si4(s3S 2 )S 32 ] s ' S2 ) (22) 



It is obvious that the matrix element of = = 1 is 

= ((Xa'XbO S ' ,S Hi|(XaXb) S ' Sz ) = 5ss'5 Sz s^s a s^5 Sb s' s (23) 

The matrix element of V^ 3 -* = Sj • Sj can be derived straightforwardly by using the recoupling 
formula Eq. fl22l . For (i, j) = (1, 2) or (3,4), the matrix element is given by 

((XA'XB'f' S '^\(XAXB) S& } = ((XA'XB'f^h ■ sMxaXb) S ' S >) 



Si 


S3 






Sl 


S3 


Sa ' 


< s 2 


Sl 


Sb' 




s 2 


Sl 


S B > 


Si 2 


S34 


S 

J 




Si 2 


S34 


S 



Sl2,S34 

-Si(Sj + 1) - Sj(Sj + 1)] 

For (2, j) = (1, 4) or (3, 2), the matrix element of V^ 3 ' = Sj • Sj is 

((xA>XB>f x \vW\(xAXB) s > s *) = ((xA'XB'f^h ■ sMxaxb) 5 * 3 *) 



(24) 
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fSs SASsSA'SB/Si 4 S32 < 


s 4 s 2 








s 4 


S2 


Sb > 


Sl4,S32 






















Sl4 S3 2 
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S14 


S32 


S 



-Sj(Sj + 1) - Sj(Sj + 1)] 
IV. Z+(4430) AND D*Di MOLECULAR STATE 



(25) 



Because the mass of Z + (4430) is close to the D*Di threshold and its width roughly is 
the same as that of Di, it is very likely that Z + (4430) is a loosely bound D*Di molecular 
state. In this section we will dynamically study whether there exists D*Di molecule state 
consistent with Z+(4430). Since m m ~ 2.422GeV, m D / ~ (2.441 ± 0.032)GeV and m D2 ~ 
2.459 GeV [361 ] . the masses of D*Di, D*D[ and D*D 2 are close to each other. Under the 
residual interaction Vj in Eq.flS]) and Eq.flHJ), these three channels would be coupled with 
each other. However, the width of D' x is very large T ~ 384 MeV 36j, consequently there 
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should be very small D*D' X component in the molecular state, otherwise it would decay so 
quickly that a weakly bound molecule can not form. As a result, we shall consider both 
D*D! and D*D 2 channels here, the effective interaction potential is induced by the pairwise 
interactions between quarks or antiquarks. Then we solve the corresponding two channels 
coupled Schrodinger equation to find whether there is bound state solutions, where we only 
concentrate on the lowest mass state. 

The model parameters employed are m u = rrid = 0.334GeV, m c = 1.776GeV, m& = 
5.102GeV, b = 0.18 GeV 2 , which is a set of fairly conventional quark model parameters. In 



Ref. 26j the screening mass /i is taken to be 0.28 GeV, which was found to be consistent with 



the string breaking mechanism and meanwhile give a good description of the charmonium 



masses 37J]. The uncertainty of screening parameter [i would be considered in the following. 
Moreover, we use a running coupling constant a s (Q 2 ), which is given by 

a » iQ2) = (33 - 2^) ln(A + QVi? 2 ) (26) 
with A = 10 and B = 0.31 GeV. Theoretical estimates for the harmonic oscillator parameter 
j3 scatter in a relative large region 0.3-0.7 GeV. Many recent quark model studies of meson 



and baryon decays use a value of (3 = 0.4 GeV 38j, |39(, therefore we assume Pa = Pb = 
0.4GeV in this work. 

In the limit that the heavy quark mass becomes infinite, the properties of the meson 
are determined by the light quark. The light quark is characterized by their total angular 
momentum, j q = s q + L, where s q is the light quark spin and L is its orbital angular 
momentum. The prime superscript (D^ or B^) is used for the state with j q = 1/2, it is very 
broad. The unprimed state (Di or B^) is used for the j q = 3/2 state, it is rather narrow. 
Heavy-light mesons are not charge conjugation eigenstates and so mixing can occur among 
the states with the same J p . The two J = 1 states Di and D' x are coherent superposition of 
the quark model 3 Pi and x Pi states 

|Di) = cos^Pi) + sin#| 3 Pi) 

|D;) = -sm^Pi) +cos#| 3 P!) (27) 

Little is known about the mixing angle 9 at present. In the heavy quark limit, the mixing 
angle is predicted to be —54.7° or 35.3° if the expectation value of the heavy quark spin-orbit 



interaction is positive or negative 



401 ] . Since the former implies that the 2 + state mass is 
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larger than that of the + state, and this agrees with the current experiment data, we shall 
employ 9 = —54.7° in the following. The above analysis applies to Bi and B' x mixing as 
well. 

There are various methods of integrating the multichannel Schrodinger equation numer- 
ically. In this work we shall employ two packages MATSCS 4l| and FESSDE2.2 42| to 
perform the numerical calculations so that the results obtained by one program can be 
checked by another. The first package is a Matlab software, and the second is written in 
Fortran 77. Both packages can fastly and accurately solve the eigenvalue problem for sys- 
tems of the coupled Schrodinger equations, and the results obtained by two codes are exactly 
the same within error. 

Calculating the relevant matrix elements of the residual interaction Vj following the 
methods outlined in section III, then we solve the coupled channel schrodinger equation 
numerically. The numerical results for the lowest energy states are listed in Table [Til We 
find that the J p = CP D*Di bound state could exist for reasonable screening mass \i. The 
binding energy is found to decrease with /i, since a smaller \i gives a stronger potential 
at short distance, which is displayed in Fig. [2j With /x = 0.28 GeV and 0.33 GeV, the 
bound state mass is about 4411.839 MeV and 4419.014 MeV respectively, and the root of 
mean square radius is 0.971 fm and 1.183 fm respectively, which are widely extended in 
space. Because the total angular momentum of S wave D*D 2 is 1, 2 or 3, it can not be 0. 
Hence D*Di will not mix with D*D2 for J p = 0~ state, then we only need to solve single 
channel Schrodinger equation in this case. Similar bound state solutions have been found 
for J p = 2~, and the binding energy is approximately the same as that of the J p = (T case 
for the same /i value. Both the 0~ and 2~ bound states are widely extended, it is a good 
feature of molecular states. The wavefunctions for the two states with \i = 0.28 GeV are 
shown in Fig. [3j For J p = 1~, bound state solutions could be found only if the screening 
mass /i is smaller than 0.16 GeV, which is quite different from the favored value 0.28 GeV. 
Therefore we tend to conclude the 1~ D*Di molecule doesn't exist. In short, both 0~ and 
2~ D*Di bound states are predicted to exist in our model, whereas only 0~ D*Di molecule 
may exist in the 7r and a exchange model from the heavy quark effective theory 16]. 

Since the production of Z + (4430) is highly suppressed in B meson decay for J p = 2~, the 
quantum J p = 0~ is favored if future experiments confirm Z + (4430) as a loosely molecular 
state. Experimentally the mass and width of Z + (4430) are fitted to be 4433±4(stat)±l(s?/st) 
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MeV and 44l 13 (stat) lf5(syst) MeV respectively. Considering the large error in the width 
measurement and the theoretical uncertainties from the screening mass /i, Z + (4430) as a 
0~ D*D! molecular state can not be excluded. More precise measurements of its mass 
and width, partial wave analysis are important to understand the nature of Z + (4430). As 
is suggested in Ref. [jjj], it is highly desirable to use the full amplitude including both 
the production and the decay processes, in performing partial wave analysis to determine 
the spin-parity of Z + (4430). If J p = _ or 2~ is favored by future partial wave analysis, 
the molecule hypothesis is strongly supported, otherwise it is not appropriate to interpret 
Z+(4430) as a D*Dx molecule. 




r(GeV ') 

FIG. 2: The potential for 0~ D*Di as a function of the separation r for different screening mass \x. 
The solid line, short dashed and dash dotted lines represent the potentials for fi = 0.23 GeV, 0.33 
GeV and 0.43 GeV respectively. 



V. APPLICATION TO OTHER HEAVY FLAVOR SYSTEMS 

In this section, we shall apply our model to 1 ++ DD* system, the heavy flavor systems 
obtained by replacing the charm quark/antiquark in Z + (4430) with bottom quark/ ant iquark 
and DD* system respectively. Possible molecular states and their static properties are stud- 
ied in detail. 
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/x(GeV) Mass(MeV) 


r rms (fm) 


P(D*D 1 ):P(D*D 2 )(%) 








44D2 438 


0.845 


100:0 


7- 


-f) 




441 1 83Q 


0.971 


100:0 






0.33 


4419.014 


1.183 


100:0 






0.16 


4427 699 


2.650 


86.747:13.253 


J= 


=1 


0.23 


no bounded 










0.28 


no bounded 












0.23 


4401.732 


0.702 


37.220:62.780 


J= 


=2 


0.28 


4414.432 


0.832 


43.148:56.852 






0.33 


4423.988 


1.272 


55.945:44.055 



TABLE II: The predictions about the mass, the root of mean square radius(rms) and the ratio of 
D*Di probability to D*D 2 probability for the bound states of the D*Di and D*D 2 system. 




(a) (b) 

FIG. 3: The radial wave functions x( r ) = rR(r) for the possible bound states of the D*Di and 
D*D 2 system, (a) and (b) respectively corresponds to J p = 0~ and J p = 2" states. 



A. 1++ DD* and X(3872) 



The narrow charmoniumlike state X(3872) was discovered by the Belle collaboration in 
the decay B + -> K+ + X(3872) followed by X(3872) -> J/^7r + 7r _ with a statistical signifi- 
cance of 10.3 a [43j]. The existence of X(3872) has been confirmed by CDF 44|, DO 45] 
and Babar collaboration 46] . the CDF collaboration measured the X(3872) mass to be 
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/i(GeV) 


Mass(MeV) r rms (fm) 


P(DD*):P(DiD 2 ) 


0.13 


3870.489 2.964 


99.966:0.034 


0.23 


no bounded 




0.28 


no bounded 





TABLE III: The predictions for the mass, the root of mean square radius(rms) and the ratio of 
DD* probability to DiD 2 probability for the bound states of the 1 ++ DD* and DiD 2 system. 

(3871.61 ± 0.16(stat) ± 0.19(sys.)) MeV. Its quantum number is strongly preferred to be 
1 ++ [471 ] . X(3872) was suggested to be DD* molecule 48], because it is very close to the 
DD* threshold. Recently several authors investigated whether a molecule corresponding to 



X(3872) could be dynamically realizable 



51 



52. 



53 



54j, notably Suzuki argued that the one 



pion exchange forces are only able to make a feeble attraction at best. We shall dynamically 
study the DD* system in our model. The 1 ++ DD* state would couple with DiD 2 under the 
pairwise residual interactions, we would like to point that the allowed quantum numbers of 
S wave D*D* are ++ , l +_ and 2 ++ . Solving the two channel coupled Schrodinger equation 
numerically, the numerical results are shown in Table IHIl We find the bound state solutions 
appear only if the screening mass \x is smaller than 0.17 GeV, with \i = 0.13 GeV the bound 
state mass is 3870.489 MeV, and it is almost completely consisted of DD*. For reasonable 
value of fi around 0.28 GeV, we can not find bound states. Consequently, X(3872) as a DD* 
molecule seems to be disfavored in our model. However, it is remarkable that unexpectedly 
large branch ratio of X(3872) — > ip(2S)^ recently was reported 55J, which indicates the mix- 
ing between the molecular state and the conventional charmonium state should be taken 
into account. This mixing effect may enhance the binding of the molecular component, this 
subject is so subtle that it is outside the scope of the present work. 



B. Bottom analog 

The bottom analog denotes the state obtained by replacing both the charm quark 
and antiquark in Z + (4430) with bottom quark and antiquark. Although the masses of P 
wave B mesons B 2 , B X) B^ and B are very close to each other j^f], |2], the widths of B' x and 
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Bo are very large. Therefore we shall only include B*Bi and B*B2 in our coupled channel 
analysis analogous to the Z + (4430) case. The matrix elements of the residual interaction Vj 
have features similar to the D*D! and D*D 2 systems except that the former is larger than 
the latter in magnitude. Furthermore, since the kinetic energy is greatly reduced compared 
with the charmed system, Z^ b should be more strongly bound than Z + (4430). The numerical 
results are shown in Table HVl It is obvious that bound state solutions with J p = 0~, \~ and 
2~ can be found, the smaller kinetic energy and deeper potential lead to two eigenstates. 
The first states are tightly bound with the binding energy from 120 to 200 GeV. The binding 
energy of the second bound states with J p = 1" and 2~ is in the range 20 to 40 GeV, and 
the B*Bx component dominates over B*B 2 . While the second state with J p = 0~ is loosely 
bound, and it disappear for \i = 0.33 GeV. The masses predicted from potential model 7| 
and QCD sum rule [10] are shown as well, our results are consistent with these predictions 
within theoretical errors. The corresponding engenstate wavefunctions with \x = 0.28 GeV 
are displayed in Fig. [5j It is remarkable that the bound state solutions predicted in our 
model are drastically different from those of the ir and a exchange model in the heavy quark 
effective theory, where only 0~ B*Bi molecule is allowed to exist 16J. Therefore experimental 
search for the bottom analog is of great interest to distinguish the different mechanisms 
in generating the molecular states. 

Clearly, the analogous state Z^ should be searched for in the Y(2S')7r + channel, where 
T(2S) can be detected by its decay into T(lS)im. Because of its large mass, at present 
the most promising place to produce Z^ b conspicuously is the large hadron colliders such as 
Tevatron and LHC. If its spin-parity is J p = l - , its neutral partner T? hh is J PC = 1 , then 
we can search for Z\ h via e + e~ annihilation at B factory. 



C. Double charged states 

The state Z^, + is obtained by replacing the charm antiquark in Z + (4430) with bottom 
antiquark, which carrying two unit electric charge. The state replacing charm quark with 
bottom quark is conjugated to Z^ + , and the static properties such as mass, rms etc are the 
same as those of Z^ + . Consequently we only need to discuss one of them, where we focus on 
Z^. + . The analysis is somewhat different from Z + (4430) and Z£ b , because the masses of D*B l5 

1(3 























/i(GeV) 


Mass(MeV) in [7] 


Mass(MeV) in [10] 


Mass(MeV) 


r rms (fm) 


P(B*B 1 ):P(B*B 2 )(%) 






0.23 






10865.468 


0.309 


100:0 












11045.211 


1.591 


100:0 


3- 


=0 


0.28 




10740 ± 120 


10886.384 


0.315 


100:0 












11049.202 


2.619 


100:0 






0.33 






10905.503 


0.322 


100:0 






0.23 






10905.552 


0.319 


38.805:61.195 












11008.389 


0.432 


61.252:38.748 


J= 


=1 


0.28 


10730 ± 100 




10922.365 


0.325 


38.431:61.569 












11013.321 


0.451 


61.637:38.364 






0.33 






10938.010 


0.333 


38.040:61.960 












11018.666 


0.477 


62.061:37.939 






0.23 






10837.374 


0.296 


24.262:75.738 












11014.714 


0.409 


76.721:23.279 


J= 


=2 


0.28 






10860.999 


0.302 


24.262:75.738 












11020.364 


0.426 


76.720:23.280 






0.33 






10882.462 


0.308 


24.276:75.724 












11026.391 


0.452 


76.828:23.172 



TABLE IV: The predictions about the mass, the root of mean square radius(rms), the probabilities 

ft n 

of the B*Bi and B*B2 components for Z^. The mass predictions in Ref. [7j and Ref. [10[] are 
shown as well. 



D*E>2, DiB* and D2B* are almost degenerate, we should solve the four channels coupled 
Schrodinger equation instead of two channels equation numerically. Since the total angular 
momentum of S wave D*B 2 and D 2 B* can not be zero, the four channels coupled Schrodinger 
equation is reduced to two channels coupled Schrodinger equations for the J p = 0~ state. 
The numerical results for the lowest states are given in Table IVl For J p = 1~ and 2~, the 
second bound state could appear for appropriate \i values. The binding energies of the first 
bound state are in the range from 70 to 100 MeV, which are larger than the binding energy 
of the D*Di system and smaller than those of Z^ b . Similar pattern has been found for the 
second bound state solution if it exists. The mass predicted in Ref. [?| from potential model 
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is shown as well, it is in agreement with our results within the theoretical errors. We plot 
the eigenstate wavefunctions in Fig. [61 This state is difficult to be produced, since both 
charm quark and bottom quark have to be produced simultaneously. The direct production 
of this state at hadron collider such as LHC and Tevatron is most promising, and we could 
search for Z^ + via the decay channel Z^, + — - > B+(2S)tt + . If the double charged state Z^, + is 
observed in future, it would be unambiguously exotic states beyond the quark model, and 
it would be a great support to the hadronic molecule picture. 







Mass(MeV) m [7| 


Mass(MeV) 


rms(im) 


P(D*B 1 ):P(D*B 2 ):P(D 1 B*):P(D 2 B*) 






0.23 




7659.800 


0.540 


81.453:0:18.547:0 


J= 




0.28 




7672.663 


0.569 


84.253:0:15.747:0 






0.33 




7683.562 


0.605 


86.883:0:13.117:0 






0.23 




7625.797 


0.447 


35.938:15.897:17.609:30.556 










7699.090 


0.640 


31.107:58.512:7.668:2.713 


J= 


=1 


0.28 


7630 ± 100 


7642.649 


0.461 


36.791:15.586:17.048:30.575 










7707.872 


0.692 


32.246:59.236:6.727:1.791 






0.33 




7657.519 


0.477 


37.799:15.270:16.504:30.426 










7715.295 


0.766 


33.693:59.341:5.924:1.042 






0.23 




7628.911 


0.450 


0.290:38.295:39.540:21.875 










7720.202 


0.921 


76.568:0.102:10.625:12.704 


J= 


--2 


0.28 




7646.478 


0.464 


0.214:37.210:40.796:21.779 










7729.487 


1.853 


88.240: 0.039: 5.549:6.173 






0.33 




7661.866 


0.480 


0.154:36.092:42.079:21.675 



TABLE V: The predictions for the mass, the root of mean square radius(rms) and the ratio between 
different components of Zi + . The mass prediction in Ref. 



7[ is also listed. 



D. DD* system 

S wave DD* system with zero isospin would be couple with D*D* under the residual in- 
teractions in Eq.([5]) and Eq.([9]), which is governed by the spin-spin interaction. In the heavy 
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quark limit, the effective potentials are induced by the interactions between two light anti- 
quarks, which is repulsive. The effective potentials for \i = 0.28 GeV are illustrated in Fig. 
HI it is obvious that the diagonal components of the effective potentials are really repulsive, 
and the off-diagonal potential is smaller than the diagonal components. Numerically solving 
the two channel coupled Schrodinger equation, we don't find bound state solutions. The 
attractive interaction is so weak that DD* bound states don't exist. The same conclusion 
has been reached from the one boson exchange model 541 ] . 




FIG. 4: The potentials for the 1=0 DD* system. The short dashed line represents the non-diagonal 
potential V\2{r), the solid and dash dotted lines respectively denote the diagonal potential Vii(r) 
and F 2 2(r). 



VI. SUMMARY AND DISCUSSIONS 

We have dynamically studied Z + (4430) and analogous heavy flavor states in quark model. 
The proximity of Z + (4430) to the D*Di threshold strongly suggests that it may be a molec- 
ular state. For a loosely bound molecule, the interaction between the constituents of the 
interacting hadrons occurs at a relatively large separation. As a consequence, the interac- 
tion will be subject to screening due to the production of dynamical quark and gluon. The 
effective charge turns out to properly describe the interactions between the constituents of 
the two hadrons at large distance, which is incorporated in this work. Because the spin- 
spin interaction is known to be important in the non- relativist ic quark models, we have 
included the spin-spin interaction in addition to the screened color-Coulomb and screened 
linear confinement interactions in our model. 

The residual interactions between two hadrons induce state mixing effect, which is taken 
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into account by solving the coupled channel Schrodinger equation numerically, where the 
second order perturbation theory can not be used anymore. We have focused on the nearly 
degenerate channels, which is a good approximation. The numerical calculations are per- 
formed with the help of MATSCS and FESSDE2.2 packages, and the results obtained by 
the two packages are the same within error. 

For the D*Di system coupled with D*D2, J p = 0" and 2~ bound states exist for reason- 
able parameter values. However, J p = 1~ bound state solution could be found only if the 
screening mass \x is smaller than 0.16 GeV. We suggest that the most favorable quantum 
is 0~, if Z + (4430) is confirmed to be a loosely bound state by future experiments. More 
precise measurements of Z + (4430) mass and width, partial wave analysis are helpful to un- 
derstand its nature. If partial wave analysis favors the quantum number - or 2~, it would 
be a strong support to the hypothesis of Z + (4430) as a molecular state. Before concluding 
that Z + (4430) is a D*Dx molecule unambiguously, we should further study the decay and 
production properties under the hadronic molecule ansatz, then compare the theoretical 
predictions with experimental data. In addition, other effects such as cusp [l^] etc should 
be taken into account, which is beyond the scope of the present work 5^|. 

The bottom analog Z^ and Z^, + are considered as well. The former is obtained by 
replacing both the charm quark and antiquark in Z + (4430) with bottom quark and antiquark, 
and the latter by replacing the charm antiquark with bottom antiquark. The second bound 
state may appear because of the smaller kinetic energy and deeper potential. The masses 
predicted in our model are in agreement with predictions from the potential model {t| and 
the QCD sum rule [lOj. We suggest to search for these states at Tevatron and LHC via 
Z+ -> T(2S)vr + and Z++ -> B+(2,S)7r + respectively. 

We have applied our model to DD* and DD* systems as well. 1 ++ DD* bound state 
solution coupled with D1D2, can be found only if the screening mass /1 is smaller than 0.17 
GeV. The mixing between the molecular state and the conventional charmonium should be 
considered to understand the nature of X(3872). For the exotic DD* system, the diagonal 
components of the effective potentials are repulsive, and the magnitude of the off-diagonal 
potentials are not attractive enough to lead to bound states. 

Our model is different from one boson exchange model and other hadronic molecule 



models 
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541 ]. consequently the predicted bound state solutions are drastically 



different from each other. We suggest that the search for the bottom analog Z^ is crucial 
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in distinguishing the different models. 

In this work, we have studied the (Qq)-(qQ) system, where Q denotes heavy quark, and q 
represents light quark. Under the short distance interactions such as the one gluon exchange 
induced constituent quark interchange interaction s |35l| . the (Qq)-(qQ) configuration may 
mix with the (QQ)-(qq) configuration. In Ref. 2l|, Swanson considered both the long 
distance one pion exchange and the short distance quark interchange interactions, He found 
the probability of the mixing of J/tpu> with DD* ranges from zero to a maximum mixing 
of 17%. For the heavy flavor systems considered in the present work, the mass difference 
between {Qq)-{qQ) and {QQ)-{qq) is larger than that in Ref. 2l|], therefore the mixing 
between the two configurations should be smaller. We expect this mixing effect plays a 
minor role here. 

The discovery of the Y(4260) and Y(4360) represents an overpopulation of the expected 



1 charmonium states. As is suggested in Ref. [50|, |57j, |58|, |59j, a possible way of recon- 
ciling Y(4260) and Y(4360) is as follows: Y(4260) is a DDi molecule, whereas Y(4360) is a 
charmonium hybrid, it is interesting to investigate whether the DDi system admits a T" 



molecular state with mass about 4260 MeV along the same line. 
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APPENDIX A: THE SPATIAL MATRIX ELEMENTS INVOLVED IN THE 
WORK 

Firstly we give the matrix elements of screened color Coulomb, screened linear confine- 
ment and spin-spin hyperfine interactions between the ground states 

Vff{ij,r) = (0,0;0,0|V; (1) (^)|0,0;0,0) = e a >^{e^[l + Erf (J- - /^)] 
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-e^[l-Erf( 



2a. 



i.i 



V$\ij, r) = (0, 0; 0, 0|K (2) (^)|0, 0; 0, 0) = e a > 2 ^-{(r - 2^a%) e~^[l + Erf (— - /i^)] 



+(r + 2/xa?,) e" r [l-Erf( 



2a 



'J 



)]} 



V$(ij,r) = (0,0;0,0|K (3) (r 4:; )|0,0;0,0) 



37TV2 



-e i j 



m i m j a-- j 



(Al) 



where a 



fl(ij) , /£(«) 



4/3- 



+ 



1/3 



5—) 



and /3s are respectively the harmonic oscillator parameters 



of the A meson and B meson. Erf(x) is the error function Erf(x) = ^= J Xe_t dt. Other 
matrix elements can be expressed in terms of the derivative of Vqq (ij, r) with respect to r, 
concretely they are given as follows 

l,l;l,-HK (fc) M|0,0;0,0> = i^^\^\ij,r) 
1, 0; 1, 0|^(r y )|0, 0; 0, 0) = J^f^ ^-vff (ij, r) 



2PaP 



B 



dr 2 



1, -1; 1, IW^M 0; 0, 0) = U l j) Jf j) -Av£\ij, r) 



0,0;l,l|K (fc) (r y )|0,0;l,l) 
1,1; 0,01^(^)10,0; 1,1) 
0,0; 1,01^(^)10,0; 1,0) 



2{3 A p B r dr 
Mij)Mij)ld V<^i j ,r) 



2p A p B r dr 



1 1 mi 92 

2j3l dr 2 



l,0;0,0|K (fc) (r,)|0,0;l,0) = U ?} S f j) ^V&\ij,r) 



0,0;l,-l|y r ( fc )(r y )|0,0;l,-l) 



2(3 A (3 B dr 2 

, fUij) 1 d 
2[3 B r dr 



V£\ij,r) 



1, -1; 0, 0\VW { r iM 0; 1, -1) = fA f/j { ^ ~V^ r) 
0,0;l,l|K (fc) (r,)|l,l;0,0) = fA ^] f f j) ~vS\ij,r) 



l,l;0,0|K {fc) (r y )|l,l;0,0) 



2(3 A (3 B r dr 

, f_MM) 1 9 

2f3\ r dr 



0,0; 1,01^(^)11,0; 0,0) = U ^J B a j^fa>) 



l,0;0,0|V;( fc) (r y )|l,0;0,0) 



2(3a(3b dr 2 
, fliij) 9 



2f3\ dr 2 



V£\ij,r) 



0, 0; 1, -1C -1; 0, 0) = UU 2 k ^ k) ~V, 



(fc) 
'oo 



(ij, r) 
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(i,-i ; o,o|K (fc) (^)|i,-i;0,o) 



2j3\ rdr- 



(A2) 
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(d) (e) (f) 

FIG. 5: The radial wave functions x( r ) = rR(r) for Z^, (a), (b) and (c) are respectively the 
wavefunctions of the first bound states with J p = 0~, J p = 1~ and J p = 2~, (d), (e) and (f) are 
the second state wavefunctions. 
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FIG. 6: The radial wave functions x{ r ) = rR{r) of Z^ + , (a), (b) and (c) are respectively the 
wavefunctions of the first bound states with J p = 0~, J p = 1~ and J p = 2~ , (d) and (e) are the 
second state wavefunctions with J p = 1~ and J p = 2~ respectively. 



27 



